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I Abstract 

, We present several families of total boolean functions which have exact quantum query com- 

plexity which is a constant multiple (between 1/2 and 2/3) of their classical query complexity, 
'(—( ', and show that optimal quantum algorithms for these functions cannot be obtained by simply 

O i' computing parities of pairs of bits. We also characterise the model of nonadaptive exact quantum 

query complexity in terms of coding theory and completely characterise the query complexity 
(Ih ' of symmetric boolean functions in this context. These results were originally inspired by nu- 

merically solving the semidefinite programs characterising quantum query complexity for small 
problem sizes. We include numerical results giving the optimal success probabilities achievable 
by quantum algorithms computing all boolean functions on up to 4 bits, and all symmetric 
boolean functions on up to 6 bits. 

> 

! 1 Introduction 

Many important quantum algorithms operate in the query complexity model. In this model, the 
quantity of interest is the number of oracle queries to the input x G {0, l}" required to compute 
some (possibly partial) function f{x). Many functions / are now known to admit quantum speed- 
up in the case where the algorithm is ahowed a constant probability of error, and the model of 
bounded-error quantum query complexity is now relatively well understood. Intriguingly, the model 
^ , of exact quantum query complexity, where the quantum algorithm must succeed with certainty, 

^ ' seems to be more mysterious. 

If / is allowed to be a partial function, it is known that there can be an exponential separation 
between exact quantum and exact classical query complexity [8], and even between exact quantum 
and bounded-error classical query complexity [4] . There are also examples of total functions where 
exact quantum algorithms allow a speed-up over classical algorithms. In particular, it is known 
that the parity of n bits can be computed exactly using only \n/2\ quantum queries [5, 10], while 
any classical algorithm (exact or randomised) must make exactly n queries. However, this is the 
best quantum speed-up known for the exact computation of total boolean functions. Even worse, 
to the authors' knowledge this is essentially the only exact quantum speed-up known. Some years 
ago, van Melkebeek, Hayes and Kutin [12] gave a quantum algorithm that computes the majority 
function on n bits exactly using only n + 1 — w{n) queries, where w{n) is the number of Is in 
the binary expansion of n, but the only quantum ingredient in this algorithm is computing the 
parity of 2 bits using 1 query. The same applies to quantum algorithms presented by Dubrovska 
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and Mischenko-Slatenkova [9], and also algorithms by Vasilieva [19, 20]. In the other direction, 
it is known that the separation between quantum and classical exact query complexity can be at 
most cubic [15] , whereas the best known relationship between bounded-error quantum and classical 
query complexity is a 6th power [3]. 

1.1 Our results 

We show that the model of exact quantum query complexity is richer than it may have hitherto 
appeared. Our results can be divided into analytical and numerical parts. On the analytical side, 
the main results are as follows. 

• We present several new families of boolean functions / for which the exact quantum query 
complexity of / is a constant multiple (between 1/2 and 2/3) of its classical query complexity, 
and we show that optimal quantum algorithms for these functions cannot be obtained by 
simply computing parities of pairs of bits. These separations are based on concatenating 
small separations obtained for functions on small numbers of bits; indeed, we give optimal 
exact quantum query algorithms for every boolean function on 3 bits. 

• We give a simple and explicit quantum algorithm for determining whether a 4-bit string has 
Hamming weight 2, using only 2 queries. Again, this cannot be done merely by computing 
parities of pairs of input bits. More generally, we give an exact algorithm which distinguishes 
between inputs of Hamming weight n/2 and Hamming weight in the set {0, l,n — l,n}, for 
all even n, using 2 queries. This generalises the well-known Deutsch-Jozsa problem [8] of 
distinguishing Hamming weight n/2 from Hamming weight in {0, n}. 

• We characterise the model of nonadaptive exact quantum query complexity in terms of a 
coding-theoretic quantity. In this setting, all queries to the input must be made up front, in 
parallel. In contrast to the classical case, there exist non-trivial quantum algorithms in this 
model. Using our characterisation result, we completely determine the nonadaptive exact 
quantum query complexity of symmetric boolean functions. 

These analytical results were inspired by numerical investigations in which we evaluated the 
best success probability achievable by quantum algorithms computing all boolean functions on 
up to 4 bits, and all symmetric boolean functions on up to 6 bits. This was achieved using the 
semidefinite programming approach of Barnum, Saks and Szegedy [2] to formulate semidefinite 
programs (SDPs) giving the precise optimal success probability for quantum algorithms using up 
to k queries, for all k < n. We then solved these SDPs numerically using the CVX package [11] for 
Matlab; the results are given in Section 6 and Appendix A. Given an SDP solution, one can then 
write down an explicit quantum query algorithm with the same parameters; we discuss how this 
can be done in Section 4. Our analytical results were based on inspecting these algorithms. 

Some further highlights from the numerical results are as follows. 

• We conjecture that the exact quantum query complexity of the two problems of determining 
whether an n-bit string has Hamming weight exactly k is precisely max{k,n — k}. We also 
conjecture that determining whether an n-bit string has Hamming weight at least k has exact 
quantum query complexity max{/c, n — A; + 1} for k > 1. Both conjectures hold numerically 
for all n < 6. 
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• In particular, we present strong numerical evidence that the algorithm of van Melkebeek, 
Hayes and Kutin [12] is not always optimal, by showing that there exists an exact quantum 
algorithm for computing the majority function on 5 bits using only 3 queries, while their 
algorithm would require 4 queries. 

• We show numerically that all boolean functions on up to 4 bits, with the exception of functions 
equivalent to the AND function, have an exact quantum query algorithm using at most 3 
queries. 

Note that Reichardt and Spalek have previously solved related SDPs (known as the adversary 
and generalised adversary bounds) for all boolean functions on up to 4 bits [18]. These SDPs 
give lower bounds on bounded-error quantum query complexity but do not characterise it exactly, 
although the generalised adversary bound does so up to a constant factor [17]. 

1.2 Organisation 

The remainder of this paper is organised as follows. After formalising some definitions in Section 
2, in Section 3 we move on to techniques for finding separations between exact quantum query 
algorithms, classical algorithms, and quantum algorithms computing parities of input bits. We 
then discuss in Section 4 how the Barnum-Saks-Szegedy SDP can be solved for small problems to 
give explicit quantum query algorithms. Section 5 gives our algorithm for for determining whether 
a 4-bit input has Hamming weight 2. In Section 6 we give optimal exact quantum query algorithms, 
found by semidefinite programming, for every boolean function on 3 bits. Our results characterising 
nonadaptive exact quantum query complexity are in Section 7, after which we conclude with a 
discussion of open problems. Two appendices detail our numerical results for all 4-bit boolean 
functions, and all symmetric functions on up to 6 bits, and also give example CVX source code. 

2 Definitions 

2.1 Generalities and boolean functions 

For any bit-string x, \x\ will denote the Hamming weight of x, and x will denote NOT(x). Cj will 
denote the bit-string with 1 in the z'th position, and elsewhere. We use the convention [X] for a 
quantity which evaluates to 1 if the statement X is true, and otherwise. We will be interested 
in boolean functions / : {0,1}" — )• {0,1}, and particularly interested in the following families of 
functions, all on n bits. PARITY is the function PARITY(j;) = [|x| is odd]. MAJ is the majority 
function where MA J (x) = > n/2]. The EXACTfc function is defined by EXACTfc(x) = [|x| = k]. 
NAE ("not- all-equal") evaluates to if all the input bits are equal, otherwise evaluates to 1. Finally, 
the threshold function Th.k{x) evaluates to 1 if and only if |x| > k. All of these are examples of 
symmetric boolean functions; a boolean function / is said to be symmetric if f{x) only depends 
on A non-symmetric function on 3 bits which we will consider is SEL, where SEL(xi, X2, xs) is 
defined to be equal to X2 if xi = 0, and equal to X3 if xi = 1. 

We say that two boolean functions / and g are isomorphic if they are equal up to negations 
and permutations of the input variables, and negation of the output variable. This relationship is 
sometimes known as NPN-equivalence. 
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2.2 Query complexity model 



An exact classical query algorithm to compute a boolean function / is given by a decision tree, or 
in other words a rooted binary tree whose vertices represent variables Xi to be queried, and whose 
edges represent branching depending on whether Xj = or = 1. The number of queries used is 
the depth of the decision tree; the minimal depth over all decision trees is the exact classical query 
complexity (aka decision tree complexity) D(f). 

We follow what is essentially the standard quantum query complexity model (see e.g. [2, 14]). 
A quantum query algorithm to compute a boolean function / : {0, l}*^ — )• {0, 1} is specified by 
a sequence of unitary operators UQ,...,Ut which do not depend on the (initially unknown) in- 
put X. These unitaries are interspersed with oracle calls Ox (which do depend on the input x). 
The final state of an algorithm that makes t queries, before the final measurement, is given by 
UfOxUt-iOx ■ ■ ■ OxUo\0) . The overall Hilbert space Ti used by the quantum query algorithm is 
split into three subspaces Tii^ "Hwork ^5 ^out- For boolean functions, "Hout is a single qubit, 
whereas the workspace Tiwork is of arbitrary size. The oracle Ox acts only on T-li^ by mapping 
\i) I— 7- (— l)^'|i) for some hidden bit string x. Tiin is n-|- 1 dimensional and indexed by basis vectors 
|0), . . . , |n); a query to xq always returns 0. The final step of the algorithm is always simply to 
measure the ^out register and return the outcome. We say that the algorithm computes / within 
error e if, on input x, the algorithm returns f{x) with probability at least 1 — e for all x. The 
exact quantum query complexity QE{f) is the minimum number of queries used by any quantum 
algorithm which computes f{x) exactly for all x. 

Note that if boolean functions / and g are isomorphic, D{f) = D(g) and Qsif) = Qsig)- 

3 Separating exact quantum and classical query complexity 

We observe that any fixed function demonstrating a separation between exact quantum and classical 
query complexity, even a small additive constant, gives rise to a constant factor multiplicative 
separation for an asymptotically growing family of functions as follows. 

Proposition 1. Let f : {0,1}'^ — t- {0,1} be a boolean function on k bits such that QeH) = "^q 
and D{f) = nic for some niq < rric- Also let g : {0, 1}*^ — t- {0, 1} be a boolean function such that 
D{g) = n. Define the family of functions fn '■ {0,1}"^ — )• {0,1} as follows: divide the nk input 
bits into blocks hi, ... ,bn of k bits each, and set fn{xi, . . . ,Xnk) = gif{bi), f{b2), . . . , f{bn))- Then 
QEifn) < ruqu and D{fn) = nicU. 

Proof. An exact quantum query algorithm for /„, using rUqU queries can be obtained simply by 
computing fibi) for each block bi individually, and then classically computing g{f{bi), . . . , f{bn)), 
so QEifn) < mqn. On the other hand, any classical algorithm for /„ must know the value of f{bi) 
for all i with certainty. Otherwise, this would imply an algorithm that computes g using fewer than 
n queries, contradicting D{g) = n. Hence D{fn) = nicn. □ 

Note that, in contrast to the classical case, it is not necessarily true that the reverse inequality 
QEifn) > mqU holds; a counterexample is provided by the PARITY function xi © X2 on 2 bits. 
Natural examples of functions g to which Proposition 1 can be applied are AND and OR. Thus 
an example of an exact query complexity separation we obtain from our results on small boolean 
functions is the following (see Section 5). 
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Theorem 2. Let EXACT^ be the boolean function on M bits defined as follows. Split the input x 
into consecutive blocks bi, . . . ,b£ containing 4 bits each, and set EXACT\{xi, . . . , X4^) = I if each 
block bi contains exactly two Is. Then QEiEXACT^) = 2£ and D{EXACT%) = AL 

Note that the lower bound (5e(EXACT2) > 21 in this theorem does not follow immediately 
from Proposition 1, but can be shown using polynomial degree arguments [3]. 

3.1 Quantum algorithms based on parity queries 

It is well-known that quantum computers can evaluate the parity of two input bits exactly using 
only one query [5] . Thus a non-trivial class of exact quantum query algorithms consists of classical 
decision trees, each of whose nodes is a query either to an individual bit of the input, or to the parity 
of two input bits. One can put a lower bound on the number of queries used by such algorithms 
(indeed, a more general class of algorithms) as follows. 

Proposition 3. Let f : {0, 1}" — )• {0, 1} be a boolean function, and let d be the degree of f as 
an n-variate polynomial over ¥2. Then any decision tree which can query the parity of any subset 
of the input variables at unit cost must make at least d queries to the input to compute f with 
certainty. 

Proof. We will show by induction that any decision tree on parities which has depth D gives rise 
to a degree D polynomial over F2. The function computed by any such tree can be written as 
pTq + (1 + p)Ti for some degree 1 polynomial p over F2 and decision trees Tq, Ti of depth at most 
D — 1. Therefore, the degree of the polynomial obtained by repeating this procedure recursively is 
at most D. If the tree computes / on every input, this polynomial must be equal to /, and hence 
be degree d. Thus D > d. □ 



4 Quantum query algorithms from semidefinite programming 

In this section we discuss, following [2], how quantum query complexity can be evaluated as a 
semidefinite programming (SDP) problem, given in Definition 1 below. In this definition, o is the 
Hadamard (entrywise) product of matrices. The following important characterisation will be the 
key to many of our results. 

Theorem 4 (Barnum, Saks and Szegedy [2]). There is a quantum query algorithm that uses t 
queries to compute a function f : {0, 1}" — {0, 1} within error e if and only if the SDP of Definition 
1 is feasible. 

Therefore, if one minimises e subject to these semidefinite constraints, one obtains the lowest 
possible error that can be achieved by a quantum algorithm which computes / using t queries. 

4.1 A prescription for quantum algorithms 

It is perhaps not immediately obvious how, given a solution to the semidefinite program of Definition 
1, to produce a quantum query algorithm with the same parameters. This was implicit in [2]; we 
now spell out explicitly how it can be done. We will use the following standard lemma from linear 
algebra (see e.g. [13]). 
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Definition 1 (Quantum query complexity primal SDP). 

For a given boolean function / : {0, 1}" — t- {0, 1} and a given integer t, find 
a sequence of 2"-dimensional real symmetric matrices {M^^^), where < i < n and 
< J < t — 1, and 2"-dimensional real symmetric matrices Fq, Ti, satisfying the 
following constraints. 



n 

n 
1=0 


— Eq 




(1) 


n 

= 

i=0 


o m/^'~^^ (for 1 < i < t - 1) 


(2) 


To + ri 


n 

= 

i=0 




(3) 



FooFo > (l-e)Fo (4) 
F.oFi > (l-e)Fi (5) 



where Ei is the matrix {x\Ei\y) = (— and Fq and Fi are diagonal matrices 
defined by {x\Fz\x) = 1 if and only if f{x) = z, and {x\Fz\x) = otherwise. 



Lemma 5. Let S = (IV'i)) ^i^-d T = be two sequences of m vectors of the same dimension. 

Define ^ = \ipi){i\, ^ = \4'j){j\- Then there is a unitary U such that U\(j)i) = iV'i) for all i if 
and only z/^^^' = $^^<I>. If such a U exists, it can be written down as follows. Let V and W be any 
isometrics satisfying ^' = VVW^, ^ = WV (i.e. isometrics occurring in polar decompositions 
of , <^), and complete V and W to unitary matrices V and W . Then U = V'{W')^. 

(i) 

Given a set of matrices M- as in Definition 1, one can derive an explicit quantum query 
algorithm completely mechanically, as follows. Use a workspace ^work of n qubits, and ignore the 
output qubit for the time being. Let the initial state be |0)|0), and let the state of the system 
at time j (i.e. after j queries have been made, and just before the {j + l)'st query is made) be 
IV'i"'^) = Y17=o where \ip'^l) is a subnormalised state in the Hilbert space Tiwork- 

(i) 

We will define the state at time < j < t — 1 in terms of the matrices M- occurring in a 
solution to the SDP by setting 

IV'S) = \fW^\x), (6) 

where \J'^vI^ is the unique positive semidefinite square root of M^'^\ It is perhaps not immedi- 
ately clear that following this prescription leads to iV'i"'^) being normalised, let alone the sequence 
IV'i'^^)) iV'i^^)) • • • ) IV'i* corresponding to a valid quantum query algorithm for all x; however, we 
will now see that this is indeed the case. 

For any 1 < j <t, define 

n 

i#) = E(-i)''i^)i^Sr'^ 

j=0 
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and set {(j)^^) = |0)|0). Also define the matrices obtained by concatenating the vectors as columns, 

xe{0,l}" a:e{0,l}" 

The vectors \(j)x^) represent the state of the system immediately after the j'th oracle call, and 
is the initial state of the system. We would like to find unitaries C/q, • • • , C/t-i mapping 
I'/'x^) I— ^ {"ipx^) for all X. By Lemma 5, such a sequence of unitaries will exist if 

($0'))t^(j) = (${i))t${i)_ 
But observe that for any < j <t — 1, 



j=0 1=0 

so (^'0))t^0) = Y17=oM^^^- Similarly, for any l<j<t, 

n n 
i=0 1=0 

= {x\(j2E,oMt'^yy). 

As, by constraint (2) in the SDP, Yli=o ^-^^ = Ya=o Ei o M^^~^^ for ah 1 < j < t - 1, this implies 
by Lemma 5 that for each j > 1 there exists a unitary Uj such that Uj\(f)x^) = \'4^x^) for all x, and 
this Uj can be determined explicitly from Lemma 5. In the case j = 0, ($(j))t$0) = Eq, and hence 

SDP constraint (1) implies the existence of a Uq such that Uq\4)^x^) = iV'i"'') for all x. 

The final constraint we need to satisfy is that the algorithm outputs the correct result. Define 
the final state of the system on input x (just before the output qubit is measured) to be 



\lx) = |0)m(\/ro|x))work|0)out + 1 0)m ( |x) ) work| l)out • 

Then 

{^x\iy) = {x\To\y) + {x\T^\y) = {x\ (r.^EioMt''^^ \y) 

by SDP constraint (3), so by a similar argument there exists a Ut such that Ut\4>x^) = |7x) for all x. 
Measuring the output qubit gives the answer with probability (x|ro|x), which by constraint (4) is 
at least 1 — e when f[x) = 0. Similarly, by constraint (5) we obtain the answer 1 with probability 
at least 1 — e when /(x) = 1. 

Observe that we have some freedom in our choice of states |^^^]); while eqn. (6) gives one 

choice which always works, it would suffice to pick any states such that {'^^i\\'4^y\) = {x\Ml^\y). In 

particular, if the rank of M- is upper bounded by r for all i, j, one can choose states of dimension 
r throughout. This would reduce the size of the Tiwork register from n qubits to [log2 r~\ qubits. 
Also observe that without loss of generality all states and unitaries occurring in a quantum query 
algorithm can be taken to be real. 
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5 EXACT2 



We now give a simple and explicit quantum algorithm for evaluating the EXACT2 function on 4 
bits using only 2 quantum queries. This algorithm was originally inspired by numerically solving 
the SDP discussed in the previous section. The algorithm docs not use any workspace (or even 
an output register) , and hence operates solely on the 5-dimensional input register indexed by basis 
states {|0), . . . , |4)}. Define a unitary matrix U by 



/O 1 1 1 1 \ 

10 1 cj a;2 

11 u'^ to 
1 a; 1 

\1 OJ 1 / 

where co = e^^^^^ is a complex cube root of 1. We begin in the state 

i=l 

and then apply Ox, then U, then Ox again. Finally, we perform the measurement consisting of a 
projection onto the state 1-0) and its orthogonal complement. If the outcome is 1^"), we output 1, 
and otherwise 0. 

The claim is that Vx := OxUOx leaves 1-0) unchanged, up to a phase factor, when x has 
Hamming weight 2, and otherwise maps |?/') into a subspace orthogonal to \^). To see that the 
claim is correct, note first that U\ip) = |0), since 1 + a; + o;^ = 0. But for x = 0000, Ox is the 
identity. Thus 

^oooolV') = |0). 

For X = 1000, Oxl^p) = 1-0) - So Vxltp) = |0) - OxU\l). But the coeflficient of |1) in U\l) is 
zero, so Ox leaves U\l) unchanged and we have 

^1000 1 V-) = \o)-u\i). 

Similar results hold for the other weight 1 strings x. 
For X = 1100, Ox\i^) = - |1) - |2), and 

Um - |1) - |2)) = |0) - i (2|0) + |1) + |2) + (a; + + |4))) 

= 1(-|1)-|2) + |3) + |4)). 

Applying Ox once more we get 

^1100 IV') = IV')- 

We get the same result for other strings of weight 2, possibly with a phase factor. For instance 

ViooilV") = ^^IV")- 

Given a string x of weight 3, we can flip all the bits and the oracle acts identically but with a 
change of sign. Thus such a string behaves like one of weight 1, and similarly a string of weight 4 
behaves like one of weight 0. Thus, for x such that \x\ 7^ 2, 14 IV') li^s in the span of |0) and U\i) for 
i = 1,2,3,4. However (V'|0) = and {tp\U\i) = i(l + w + w^) = for i = 1,2,3,4. So this subspace 
is orthogonal to |V'), proving our claim. 
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5.1 Distinguishing weights and 1 from balanced strings 



We would ideally like to understand the number of queries required to solve the EXACT^ function 
on n input bits, for all n and k. As an intermediate goal generalising our solution for n = 4, one 
can ask for a two-query algorithm, for any even re, that distinguishes strings of weight and 1 from 
strings of weight n/2. We take the previous space, with basis vectors z = 0, 1, . . . , n, and tensor 
it with an ancilla space. In the ancilla space we select vectors |aj), i = 1, . . . , re, of unit length with 
inner products {ai\aj) = c, for i ^ j, where c is arbitrary such that |c| < 1. Such vectors can always 
be found. We also select some vector |0) in the ancilla space. 

The oracle acts on the first space, so Ox\i)\j) = (— for any Recall that is the 
string with a 1 at position i and O's elsewhere, and let b be the "balanced" string with I's in the 
first n/2 places. The algorithm starts with the state {(p) = Yl^=i 10 10) ('we keep unnormalised 
for simplicity). As before, Vx := OxUOx, with U to be defined shortly. The aim is to show that 
Vh\<j)) is orthogonal to the subspace spanned by V^Ji/*) for all i, which thus discriminates between 
balanced and weight 1 strings (discriminating weight strings will follow automatically from this) . 

Now define U by its action on the states \Ti) := Oe^\4>) = — 2\i)\0): 

U\t,) = am +13 (j2 - IOI«2^-i)j (7) 

where the subscript of the a's is taken mod re. This will be an isometry (which can be extended to 
a unitary on the whole tensor product space) if 

a2 + (re-l)/32 = re, (8) 
a2 + (re-2)/32c = re-4. (9) 

Note that U\Ti) = Ve^lcj)), since U\Ti) is invariant under Oe^- 
U can be defined on Oh\<p) by using 

OM) = 




giving, up to an overall factor of /3, 



/n/2 n \ n n n 

Vb\<l>) = -[Y.O,U\n) - O.C/h) = j;j;7r,fc|j)|afc+,_i) + ^|i)|a2,- 

\i=l i=n/2+l j j=lk=l j=l 



(10) 



where 



TTjk = -1 if i e L,j e L, TTjk = 1 if i e L,j e H, 

TTjk = I if i £ H,j £ L, TTjk = -I if i £ H,j £ H, 

and L = [l,n/2],H = [re/2 + l,re]. Combining (10) with (7) gives 
{<P\VbVei\(l)) = 'Y'^-n-jkia-j+k-ilaj+i-i) - ^(ai+j-i|«2i-i) + ^(a2i-i|ai+j_i) - (a2i_i|a2i-i) 

j k j j 

= 0-(c-l) + (l + (re-l)c)-l 
= 1 + (re - 2)c. 
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Thus for orthogonality we require n = — l/(n — 2), which, with (8) and (9), gives = (n - 2)2/n 
and = 4/n. 

FinaUy, we would like to show orthogonality for the weight zero sequence, i.e. orthogonality of 
Vb\<j)) and U\(l)). However, this follows from what we proved above, using 

1 " 

2 = 1 

together with the fact that U\Ti) is invariant under Oe^- 

6 Exact quantum query algorithms for small functions 

Having whetted our appetite with the EXACT2 problem, we now turn to quantum algorithms for 
other small boolean functions. For each function on n bits we considered, we calculated, via the 
SDP of Definition 1, the best possible success probability achievable by quantum algorithms making 
t queries, for t = 1, . . . , n — 1 (any function can clearly be computed exactly using n queries). We 
did this for all functions on up to 4 input bits, and for all symmetric functions on up to 6 bits. We 
used the convex optimisation package CVX [11] for Matlab, which allows optimisation problems 
to be specified in a simple and intuitive way; see Appendix B for source code. The CVX package 
allows the choice of underlying solvers SeDuMi and SDPT3. We used SeDuMi for the results given 
below, and also checked the results with SDPT3, which gave the same values up to a difference of 
at most 0.001. The numerical results for functions on 4 bits, and symmetric functions on up to 6 
bits, are deferred to Appendix A. 

Note that there is a basic issue with calculating exact quantum query complexity numerically, 
which is that one receives a numerical solution from the SDP solver, which is not exact. If the SDP 
solver claims that there exists a quantum query algorithm that computes some function / using k 
queries with success probability at least 0.999, for example, one cannot be sure that this algorithm 
is actually exact. In the case of all functions on up to 3 bits, we therefore give explicit optimal exact 
quantum query algorithms. These algorithms were obtained by a somewhat laborious process of 
taking the numerically obtained (real- valued, approximate) solutions to the SDP and using these 
as a guide to find exact solutions. 

For completeness, we begin by giving optimal exact quantum query algorithms for all boolean 
functions of 1 and 2 bits. In what follows, the tables are indexed by function ID; the ID of each 
function is the integer obtained by converting its truth table from binary. Columns give the optimal 
success probability that can be achieved by quantum algorithms making 1, . . . , n — 1 queries. Entries 
are starred when there is a nonadaptive exact quantum algorithm using that number of queries (see 
Section 7). 

6.1 Functions of up to 2 bits 

Up to isomorphism, the only non-constant function of 1 bit is f{xi) = xi, which clearly requires 
exactly one query. In the case of 2 bits, there are two classes of functions. 



ID 


Function 


1 query 


1 


Xi A X2 


0.900 


6 


Xl © X2 


1* 
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An optimal quantum algorithm for the function xi (Bx2 proceeds as follows [5]. Input the state 
-^(|1) + |2)) into the oracle to produce + {—lY'-^\x2))- Perform a Hadamard gate 

(with respect to the basis {|1), |2)}), measure in the basis {|1), |2)}, and output if "1" is measured, 
and 1 if "2" is measured. It is easy to see that this algorithm succeeds with certainty. 

6.2 Functions of 3 bits 

The following table lists all boolean functions depending on 3 bits, up to isomorphism. 



ID 


Function 


1 query 


2 queries 


F2 degree 


D(f) 


1 


xi A X2 A X3 


0.800 


0.980 


3 


3 


6 


Xi A (X2 © X3) 


0.667 


1* 


2 


3 


7 


xi A {x2 V X3) 


0.773 


1 


3 


3 


22 


EXACT2 


0.571 


1 


3 


3 


23 


MAJ 


0.667 


1 


2 


3 


30 


Xi © (X2 V X3) 


0.667 


1 


2 


3 


53 


SEL(xi,X2,X3) 


0.854 


1 


2 


2 


67 


(xi A X2) V (xi A X2 A X3) 


0.773 


1 


3 


3 


105 


PARITY 


0.500 


1* 


1 


3 


126 


NAE 


0.900 


1* 


2 


3 



Observe that the AND function requires 3 queries to be computed exactly; in fact, it has been 
known for some time that AND on n bits has Qe(AND) = n [3]. For most of the other functions, 
an optimal exact quantum algorithm is easy to determine, based only on classical queries and 
computing the parity of two bits using one query. 

• xi A (x2 © X3): Query xi and evaluate X2 © X3 using one query. 

• MAJ: First evaluate xi © X2. If the answer is 1, then output X3, otherwise output xi. This 
works because if xi and X2 are different, then there will be at least two I's in total if and only 
if X3 is 1. If xi and X2 are the same, there will be at least two I's if and only if xi is 1. 

• xi © (x2 V X3): This function is equivalent to (X3 A X2) V (X3 A (xi © X2)). So query X3 first, 
then either query X2 or xi © X2. 

• SEL(xi, X2, X3): Query xi first, then either output X2 or X3. 

• PARITY: Evaluate xi © X2, query X3, and take the exclusive OR of the two. 

• NAE: This function is equivalent to (xi © X2) V (xi © X3). 

However, the three remaining functions (EXACT2, xi A (x2 V X3) and (xi A X2) V (xi A X2 A X3)) 
do not have such straightforward optimal algorithms. Indeed, by Proposition 3, they cannot be 
computed using 2 queries by any algorithm which is a decision tree on parity queries. 

In the case of EXACT2, we obtain an optimal algorithm by appending an additional zero bit 
and computing EXACT2 on 4 bits (see Section 5). We now give quantum query algorithms for the 
two remaining functions. Rather than writing out the unitary operators arising in the algorithm 
explicitly, we simply give expressions for matrices forming a exact solution to the query complexity 
SDP. We stress that, given these matrices, one can follow the procedure of Section 4 to find an 
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explicit quantum algorithm completely mechanically. The matrices are fully specified by their 
non-zero eigenvalues and eigenvectors; note that, for readability, we have neither normalised nor 
orthogonalised the eigenvectors. 

6.2.1 Xi A (X2 V .Tg) 



A/r 'i' T* 1 "V" 




Th',! o*OTi"\7or'"i"r~»T'c: 

J_Jl&dl V LUl o 


A^(0) ,^(0) ,^(0) ,,(0) 

Mq ',Ml ',M^ ',M^ 


2 


(1,1,1,1,1,1,1,1) 




3/2 


(1. 1, 1. 1,0,0,0,0) 




1 


{(-1, 1,-1,1,0,2,0,2), (0,-1, 1,0,0,-1, 1,0)} 




1 


(1,0,0,-1,2,1,1,0) 




1/2 


(-1,0,-1,0,0,1,0,1) 




1 


(1,0,0,-1,2,1,1,0) 




1/2 


(-1,-1,0,0,0,0,1,1) 




3/4 


{(0,1, 0,1, 0,1, 0,1), (0,-1, 1,0, 0,-1, 1,0)} 


To 


5/2 


(3,2,2,3,2,0,0,0) 




1 


{(0, -1, 0, 0, 1, 0, 0, 0), (0, -1, 1, 0, 0, 0, 0, 0)} 




1/2 


(-1,0,0,1,0,0,0,0) 




3/2 


{(0,0, 0,0, 0,1,0,1), (0,0, 0,0, 0,-1,1,0)} 



6.2.2 (xi A X2) V [xi A X2 A X3) 



Matrix 


Eigenvalues 


Eigenvectors 




2 


(1,1,1,1,1,1,1,1) 




1 

3/4 
1/4 


(-2,-1,-1,0,-1,0,0,1) 
(0,0,-1,-1,2,2,1,1) 
(0,0,1,1,0,0,1,1) 




1 

3/4 
1/4 


(-2,-1,-1,0,-1,0,0,1) 
(0,-1,2,1,-1,-2,1,0) 
(0,-1,0,-1,1,0,1,0) 




1 

3/4 
1/4 


(0,1,1,2,1,2,2,3) 

(0,-1,0,-1,1,0,1,0) 
(0,-1,2,1,-1,-2,1,0) 




1 

3/4 
1/4 


(0,3,-1,2,-1,2,-2,1) 

(0,0,1,1,0,0,1,1) 
(0,0,-1,-1,2,2,1,1) 


To 


i(5 + ^/5) 
3/2 
1 

i(5-v^) 


(1 + V5,0, -1 + i(5 + V5), 1, -1 + i(5 + 1, 0,0) 

(0,0,0,-1,0,1,0,0) 

(0,0,-1,0,1,0,0,0) 
(1 - 0, -1 + ^(5 - v/5), 1, -1 + i(5 - v/5), 1, 0, 0) 




3/2 


{(0,-1, 0,0, 0,0, 0,1), (0,1, 0,0, 0,0, 1,0)} 
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7 Nonadaptive exact quantum query complexity 



We now turn to a very restricted model of exact query complexity, in which the algorithm's queries 
are required to be nonadaptive. In other words, the choice of which input variables to query cannot 
depend on the result of previous queries, so the algorithm must choose which variables to query at 
the start. For a boolean function /, let Q^{f) be the nonadaptive quantum and classical 

exact query complexities of /, i.e. the minimum number of nonadaptive queries required to compute 
/ with certainty. This model is extremely restricted classically, as we see from the following easy 
proposition. 

Proposition 6. For any boolean function f depending on n variables, D^"'{f) = n. 

Proof. A nonadaptive exact classical query algorithm A making k queries is specified by a list of 
k fixed variables which are queried. If /c < n, there must exist a variable i which is not queried, 
but on which / depends. Thus there must exist an input x such that if bit i is flipped, A does not 
notice the difference, so A cannot be correct on every input. □ 

In the case of nonadaptive quantum query complexity, things are more interesting. An optimal 
quantum algorithm for PARITY is nonadaptive and uses \n/2\ quantum queries [5, 10], so quantum 
algorithms can indeed achieve an advantage over classical algorithms in this model. However, it is 
also known that this separation by a factor of 2 is optimal for total functions in the nonadaptive 
model [16]. 

We now introduce some additional notation. For any boolean function / : {0, 1}" — t- {0, 1}, 
define 

Sf := {z : \/x, f{x) = f{x + z)}, 

where addition is over the group Zg; i.e. Sj is the set of translations of the input under which 
/ is invariant. Note that Sf is a subspace of {0,1}". For any subspace S C {0,1}", let S-^ 
denote the orthogonal subspace to S, i.e. S-^ = {x : x • s = 0,Vs G S}. Finally, let d{x,S) 
denote the maximum Hamming distance between a bit-string x G {0, 1}" and a subset S C {0, 1}"": 
d{x, S) = maxyg5 d{x, y). Then we have the following theorem. 

Theorem 7. For any boolean function f : {0, 1}" — t- {0, 1}, 

QW'if) = maxd(3;,y)= min d(x,Sf). 

xe{o,i}" yf^sf xe{o,i}" 

We have thus completely characterised the nonadaptive quantum query complexity of /. In 
the coding theory literature, the quantity minj-gjo,!}" d{x,S-f) is known as the radius of the code 
Sj- [6]. Observe that Theorem 7 implies that Q^{f) can be computed exactly in time polynomial 
in 2". We now prove this theorem and then draw some corollaries. In the proof it will be convenient 
to use the notation 

/» = ^ E (-ir'/(y) 

ye{0,l}" 

for the Fourier transform (over Z2 ) of some function / : {0, 1}" — )• M. 

Proof of Theorem 7. Any nonadaptive quantum algorithm making k queries to the input x corre- 
sponds to a choice of a state {ip), which is input to k copies of the oracle (i.e. the unitary operator 
O®^), followed by a projective measurement to determine whether f{x) = or f{x) = 1. Label 
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computational basis states by a length k string of integers {ii, . . . , ik) in the range {0, . . . , n}; each 
such string represents a list of variables queried, with representing a "null query" which does 
nothing. 0®^ acts on these basis states by mapping 

\ii,...,ik) ^ 

Now note that we can restrict ourselves to query strings in non-decreasing order, and containing at 
most one of each integer between 1 and n. The first of these is because querying any permutation 
of a string is equivalent to querying the string itself. The second is because querying the same 
index twice does nothing, and hence is equivalent to the null query. 

These strings are now in obvious one-to-one correspondence with the set of n-bit strings of 
Hamming weight at most k. Thus the state we obtain from applying Of^ to an arbitrary input 
state is of the form 

s6{0,l}»,|s|<fe 

A nonadaptive quantum query algorithm computing / exactly using k queries exists if and only if 
there exists a set {ccs} such that (V'xl'^j/) = for all x, y such that f{x) ^ f{y), or in other words 

se{o,i}",|s|<fe 

For brevity, write w{s) = la^p. The above constraint says that, for all z ^ Sj, 

se{0,l}",|s|<A: 

Considering w{s) as a function w : {0, 1}" — >■ M such that w(s) = for \s\ > k, in Fourier-analytic 
terminology the constraint says that 

w{z) = Oifz^Sf, 

i.e. that w is only supported on the subspace Sf. This is equivalent to the constraint that 

w{s) = - — -j— w{s + t) for all s. 

To see this, define the function Psf{x) = [x G Sf], and note that Psf{t) = j^I* ^ S^]. Letting * 
denote convolution over Z2 (i.e. (/ * g){x) = f{y)g{x + y)), by Fourier duality we have 

PSfiji^ = w Psf * w = w 4^ ^ w{s + t)Psf {t) = w{s) for all s 

t 

<^ T^lT X] '"^(^ + *) = '"^(^) fo^ 
'•^Z ' tesf 

Thus w is uniform on cosets of S-^. As w is not identically zero, there must be a coset t + Sj- such 
that every element s e t + Sj- has Hamming weight at most k. If Sf = 0, then Sj- = {0, 1}" and 
hence has only one coset, which contains 1". Hence we must have k = n. More generally, we have 
that Q^if ) is the minimal k such that there exists a t satisfying |,s| < k for all s £ t + S'f. In other 
words, Q^{f) is the minimal k such that there exists a t satisfying d{s,t) < k for all s G S^. □ 
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We observe from this proof that it is without loss of generahty that any nonadaptive exact 
quantum algorithm can be described as picking a coset of Sj' and querying everything in that 
subset uniformly. Explicitly, we have the following algorithm. 

1. Let t e {0, !}"■ be a bit-string such that d{t, S'j ) = k. 

2. Produce the state of n qubits i^^^^^a Z^sGt+Sj^ ^ ^ queries to the oracle. 

3. Perform Hadamards on every qubit of the resulting state and measure to get outcome x. 

4. Output f{x). 

One can easily verify that in fact /(x) = /(x) with certainty. We note that this is reminiscent of 
an algorithm of van Dam [7] which learns x itself with bounded error using n/2 + 0{^/n) queries to 
the oracle. Here, we also compute a partial Fourier transform from which f{x) can be determined, 
but our algorithm succeeds with certainty. 

We now draw some corollaries from Theorem 7. 

Corollary 8. For any boolean function f : {0, 1}" — )• {0, 1} such that f is not invariant under any 
translation, Q'^{f) = n. 

Proof. If / is not invariant under any translation, Sj = $ and hence Sj- = {0,1}". For any 
bit-string x G {0, 1}", there exists a y G {0, 1}" such that d{x,y) = n. Hence Q^{f) = n. □ 

One could also have observed this corollary by noting that Q^{f) depends only on Sj, and for 
the AND function (which has (5^(AND) = n [3]), 5and = 0- The corollary implies that only an 
exponentially small fraction of boolean functions / have Q^{f) < n. One class of functions that 
do satisfy this is given by the following corollary. 

Corollary 9. For any boolean function f : {0, 1}"" — )• {0, 1} such that f{x) = f{x) for all x, 
Ql^{f)<n-l. 

Proof. The constraint f{x) = f{x) for all x is equivalent to 1") = f{x) for all x, so {0", 1"} C 

Sf, implying S'j C {x : |x| even}. If n is odd, then d{G^,S^) <n — \ (as 1" has odd Hamming 
weight, there is no even weight bit string distance n from 0"). Similarly, if n is even, ^(10""^, S'j ) < 
n- 1. Hence Q^"(/) < n- 1. □ 

An explicit nonadaptive quantum algorithm achieving this query complexity proceeds as follows. 
Evaluate 7/^ := xi©Xfc, for 2 < /c < n, at a cost of n — 1 queries in total, then output f{0,y2,---, yn)- 
If xi = 0, this is simply /(x), while if xi = 1, this is /(x) = /(x). 

Corollary 10. For any boolean function f : {0, 1}"" — ?■ {0, 1} that depends on all n input bits, 
Ql^{f)>\n/2\. 

Proof. We need to show that, for any bit-string x, we can find an element y S Sj' such that 
d{x,y) > n/2. As / depends on all its input bits, for all i £ {1, . . . ,n}, ^ Sj. This implies that, 
for all i G {1, . . . ,n}, there is at least one element of Sj' whose i'th bit is 1. Thus, if we pick an 
element y G Sj' at random, each bit of y will be or 1 with equal probability, so the expectation 
of d{x,y) is exactly n/2. Therefore, d(x,5^) > n/2. □ 
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Corollary 10 was previously proven in [16] via a different method. We can also show that 
functions whose nonadaptive exact quantum query complexity is minimal are of very restricted 
form. 

Corollary 11. Let f : {0, 1}" — )• {0, 1} be a boolean function that depends on all n input bits and 
such that Q^{f) = n/2. Then there exists an x such that d{x,y) = n/2 for all y G . 

Proof. By the proof of Corollary 10, ¥.yi^gjd{x,y) = n/2. Thus, if d{x,y) < n/2 for all y G 5^ , we 
must have d{x,y) = n/2 for all y G S'^. □ 



7.1 Symmetric boolean functions 

It turns out that we can apply Theorem 7 to completely characterise the nonadaptive exact quantum 
query complexity of symmetric boolean functions, via the following quadrichotomy. 

Theorem 12. Let f : {0, 1}" — )• {0, 1} be symmetric. Then exactly one of the following four 
possibilities is true. 

L f is constant and Q^{f) = 0. 

2. f is the PARITY function or its negation and Q^{f) = \n/2\. 

3. f satisfies f{x) = f{x) (but is not constant, the PARITY function or its negation) and 
Q%%f) = n-l. 

4- f is none of the above and Q^{f) = n. 

We will prove Theorem 12 using the following lemma, whose proof is given afterwards. 

Lemma 13. Let f : {0, 1}'" — {0, 1} be symmetric and satisfy f{x) = f{x + a) for all x G {0, 1}*^, 
for some a with 1 < \a\ < n — 1. Then, if \a\ is odd, f is constant. If \a\ is even, f is constant, 
PARITY or its negation. 

Proof of Theorem 12. First assume there is an a with 1 < |a| < n — 1 such that f{x) = f{x + a) for 
all X G {0, 1}". Then, by Lemma 13, / is constant, PARITY or its negation. If / is constant then 
clearly Q^f (/) = 0. If / is PARITY or its negation, by the result [10] of Farhi et al. Q'^Hf) = [n/2] . 
On the other hand, if there is no a with 1 < |a| < n— 1 such that f{x) = f{x + a) for all x G {0, l}", 
but there is such an a with \a\ = n, f{x) = f{x) and by Corollary 9 Q^{f) = n — 1. Finally, if 
there is no a / 0" such that f {x) = f{x + a) for ah x G {0, 1}", by Corollary 8 Qf{f) = n. □ 

It will be convenient to prove Lemma 13 using Fourier analysis, based on the following well- 
known fact. 

Fact 14. For any f : {0, 1}" M and for any a G {0, 1}", if f{x) = f{x + a) for all x G {0, 1}", 
then for all b such that |a A 6| is odd, f{b) = 0. 

Proof of Lemma 13. Note that, because / is symmetric, if f{s) = for some s with \s\ = k, f{t) = 
for all t with \t\ = k. Without loss of generality, assume a consists of j ones followed hj n — j zeroes 
(i.e. is of the form 1 ... 10 ... 0). Consider the bit-string s which is 1 on the set {1, . . . , 2A: + 1} for 
some < k < {j — l)/2, and the set {n — i + 1, . . . ,i}, for some 1 < i < j. By Fact 14, for any 
such bit-string s, f{s) = 0. By varying k and i, we can vary |s| arbitrarily between 1 and either n 
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(if \a\ is odd), or n — 1 (if \a\ is even). Thus, if \a\ is odd, f{s) = for all s 7^ 0", so / is constant. 
Otherwise, if \a\ is even, f{s) = for all s ^ {O", 1"}. The only boolean functions satisfying this 
are constant functions, PARITY and its negation. □ 

8 Open problems 

It is a very tempting conjecture that QE(EXACTj!c) = max{A;,n — k}. It is easy to see that the 
lower bound Q£;(EXACTfc) > max{A;,?i — k} holds; by setting min{A;,n — k} input bits to we 
obtain a function equivalent to the AND function on i := max{A;,n — k} bits, which has exact 
quantum query complexity i. So it would suffice to prove the upper bound (5£;(EXACT„/2) ^ ^^/^ 
for all even n to prove this conjecture. To see this, note that for any i, an algorithm for the function 
EXACT^ on n bits gives an algorithm for EXACT^ on n' > n bits, simply by appending extra zero 
bits. 

More generally, the problem of finding any boolean function / such that Q_e(/) < D{f)/2 still 
remains. We are hopeful that the numerical techniques used in this paper may prove helpful in 
resolving this question. 
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A Numerical results for functions on up to 6 bits 

In this appendix we collate our numerical results concerning the optimal success probability achiev- 
able by quantum algorithms for all boolean functions on 4 input bits, and symmetric boolean func- 
tions on 5 and 6 input bits. We split the results into sections according to the number of bits 
on which the functions depend. Note that each section on functions of k bits does not include 
functions which only depend on fewer than k bits. In the following tables, entries are starred when 
there is a nonadaptive exact quantum algorithm using that number of queries (see Section 7). An 
entry "1" means that the SDP solver claims a solution with success probability greater than 0.999; 
note that this does not strictly speaking imply the existence of an exact algorithm using that num- 
ber of queries. We use the notation SYM(co, . . . , c„) to mean the symmetric function / such that 
f{x) = c\x\- III the tables of symmetric functions, we simply identify each function with the vector 
(co, . . . ,c„). 

Note that for all non-constant symmetric /, the decision tree complexity D{f) = n, but this is 
not the case for QeH)- For most functions / on 4 bits, D{f) is easily verified to be 4 via polynomial 
degree arguments; we calculated D{f) for the remaining functions / in an ad hoc fashion. 
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A.l Functions of 4 bits 



ID 


Function 


1 query 


2 queries 


3 queries 


D(f) 


1 


xi A X2 A a;3 A X4 


0.735 


0.962 


0.996 


4 


6 




0.654 


0.931 


1* 


4 


7 




0.750 


0.954 


1 


4 


22 




0.572 


0.906 


1 


4 


23 




0.667 


0.926 


1 


4 


24 


Xi A -iNAE(x2, ^3, 2^4) 


0.654 


0.931 


1* 


4 


25 




0.640 


0.961 


1 


4 


27 


xi A SEL(a;4, X2,X3) 


0.667 


0.965 


1 


3 


30 




0.600 


0.956 


1 


4 


31 




0.718 


0.970 


1 


4 


61 




0.643 


0.976 


1 


4 


105 




0.500 


0.900 


1* 


4 


107 




0.571 


0.941 


1 


4 


111 




0.662 


0.968 


1* 


4 


126 


xi A NAE(x2, 0:3, 0:4) 


0.667 


0.947 


1* 


4 


127 




0.727 


0.972 


1 


4 


278 


EXACT3 


0.529 


0.884 


1 


4 


279 


Th3 


0.643 


0.900 


1 


4 


280 




0.572 


0.906 


1 


4 


281 




0.600 


0.956 


1 


4 


282 




0.571 


0.936 


1 


4 


283 




0.637 


0.959 


1 


4 


286 




0.546 


0.932 


1 


4 


287 




0.659 


0.945 


1 


4 


300 




0.571 


0.936 


1 


4 


301 




0.572 


0.964 


1 


4 


303 


SEL(a;3, xi A X2, SEL(x4, xi,X2)) 


0.644 


0.966 


1 


3 


316 




0.562 


0.962 


1 


4 


317 


SEL(x3, xi A X2, SEL(x2, xi,Xi)) 


0.572 


0.980 


1 


3 


318 




0.546 


0.956 


1 


4 


319 




0.640 


0.972 


1 


4 


360 




0.529 


0.884 


1 


4 


361 




0.500 


0.916 


1 


4 


362 




0.546 


0.932 


1 


4 


363 




0.546 


0.955 


1 


4 


366 




0.546 


0.956 


1 


4 


367 




0.571 


0.969 


1 


4 


382 




0.546 


0.923 


1 


4 


383 




0.600 


0.946 


1 


4 


384 


NAE(xi,a;2,a;3,iC4) 


0.800 


0.980 


1* 


4 


385 




0.750 


0.954 


1 


4 


386 




0.640 


0.961 


1 


4 


387 




0.667 


0.965 


1 


4 


390 




0.571 


0.936 


1 


4 


391 




0.637 


0.959 


1 


4 
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ID 


Function 


1 query 


2 queries 


3 queries 


D(f) 


393 


SEL(a;3, xi A X4, X2 A X4) 


0.667 


0.965 


1 


3 


395 




0.724 


0.963 


1 


4 


399 




0.751 


0.980 


1 


4 


406 




0.500 


0.916 


1 


4 


407 




0.572 


0.940 


1 


4 


408 




0.600 


0.956 


1 


4 


409 




0.643 


0.976 


1 


4 


410 




0.572 


0.964 


1 


4 


411 




0.656 


0.969 


1 


4 


414 




0.546 


0.955 


1 


4 


415 




0.642 


0.965 


1 


4 


424 




0.667 


0.926 


1 


4 


425 




0.637 


0.959 


1 


4 


426 




0.718 


0.970 


1 


4 


427 


SEL(x3, xi A X4, SEL(a;4, xi , X2)) 


0.751 


0.980 


1 


3 


428 




0.637 


0.959 


1 


4 


429 


SEL(a;2, xi A X4, SEL(x3, xi , X4)) 


0.656 


0.969 


1 


3 


430 




0.644 


0.966 


1 


4 


431 




0.710 


0.977 


1 


4 


444 




0.572 


0.980 


1 


4 


445 




0.641 


0.965 


1 


4 


446 




0.572 


0.969 


1 


4 


447 




0.667 


0.980 


1 


4 


488 




0.643 


0.900 


1 


4 


489 




0.572 


0.940 


1 


4 


490 




0.659 


0.945 


1 


4 


491 




0.642 


0.965 


1 


4 


494 




0.640 


0.972 


1 


4 


495 


SEL(x3, xi, SEL(x4, xi, X2}) 


0.667 


0.980 


1 


3 


510 




0.600 


0.946 


1 


4 


829 




0.563 


0.975 


1 


4 


854 




0.598 


0.955 


1 


4 


855 




0.714 


0.969 


1 


4 


856 




0.572 


0.964 


1 


4 


857 




0.579 


0.961 


1 


4 


858 


SEL(xi, X2 A x^, x2 X4) 


0.572 


0.980 


1 


3 


859 




0.628 


0.974 


1 


4 


862 




0.572 


0.966 


1 


4 


863 


SEL(xi, X2 A X3, x2 V X4) 


0.667 


0.986 


1 


3 


872 




0.546 


0.932 


1 


4 


873 




0.500 


0.946 


1 


4 


874 




0.598 


0.955 


1 


4 


875 




0.572 


0.951 


1 


4 


876 




0.546 


0.956 


1 


4 


877 




0.545 


0.961 


1 


4 


878 




0.572 


0.966 


1 


4 


879 




0.600 


0.966 


1 


4 


892 




0.563 


0.975 


1 


4 
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ID 


Function 


1 query 


2 queries 


3 queries 


D(f) 


893 




0.571 


0.966 


1 


4 


894 




0.572 


0.947 


1 


4 


961 




0.718 


0.970 


1 


4 


965 


SEL(x2, xi A X3, SEL(aji , X3, X4)) 


0.751 


0.980 


1 


3 


966 




0.644 


0.966 


1 


4 


967 




0.710 


0.977 


1 


4 


980 




0.659 


0.945 


1 


4 


981 




0.714 


0.969 


1 


4 


982 




0.572 


0.951 


1 


4 


983 




0.661 


0.965 


1 


4 


984 


SELfxi xoAx-i SELixA x-? xo)) 


0.644 


0.966 


1 


3 


985 




0.628 


0.974 


1 


4 


987 


SELfx4 SELfxq xi x-A SELfx9 xi x-\]\ 


0.661 


0.965 


1 


3 


988 




0.640 


0.972 


1 


4 


989 


SEL(xi, X2 A X3, V X4) 


0.667 


0.986 


1 


3 


990 




0.600 


0.966 


1 


4 


1632 




0.667 


1* 


1* 


4 


1633 




0.562 


0.962 


1 


4 


1634 




0.643 


0.976 


1 


4 


1635 




0.572 


0.980 


1 


4 


1638 




0.667 


0.947 


1* 


4 


1639 




0.641 


0.965 


1 


4 


1641 




0.500 


0.936 


1* 


4 


1643 




0.561 


0.966 


1 


4 


1647 


MAJ(xi, X2, X3 © X4) 


0.667 


1 


1* 


4 


1650 




0.656 


0.969 


1 


4 


1651 




0.628 


0.974 


1 


4 


1654 




0.641 


0.965 


1 


4 


1656 




0.546 


0.956 


1 


4 


1657 




0.500 


0.964 


1 


4 


1658 




0.571 


0.966 


1 


4 


1659 




0.571 


0.962 


1 


4 


1662 




0.600 


0.954 


1 


4 


1680 


(xi X2) A (xi © X3 © X4) 


0.500 


0.900 


1* 


4 


1681 




0.500 


0.916 


1 


4 


1683 




0.500 


0.946 


1 


4 


1686 




0.500 


0.936 


1* 


4 


1687 




0.500 


0.964 


1 


4 


1695 


SEL(x3 © X4, xi , X2) 


0.500 


1 


1* 


3 


1712 




0.571 


0.941 


1 


4 


1713 




0.546 


0.955 


1 


4 


1714 




0.572 


0.940 


1 


4 


1715 




0.572 


0.951 


1 


4 


1716 




0.546 


0.955 


1 


4 


1717 




0.545 


0.961 


1 


4 


1718 




0.561 


0.966 


1 


4 


1719 


SEL(x4, SEL(x2, xi, X3), SEL(x3, X2, xi)) 


0.572 


0.962 


1 


3 


1721 




0.500 


0.964 


1 


4 
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ID 


Function 


1 query 


2 queries 


3 queries 


D(f) 


1725 




0.529 


0.955 


1 


4 


1776 




0.662 


0.967 


1* 


4 


1777 




0.572 


0.969 


1 


4 


1778 




0.642 


0.965 


1 


4 


1782 


SEL(x2, Xi, X3 © X4) 


0.667 


1 


1* 


3 


1785 


x-i ffi (x9 A (x-i ffi Xd)) 


0.500 


1 


1* 


4 


1910 




0.600 


0.954 


1 


4 


1912 




0.546 


0.923 


1 


4 


1913 




0.529 


0.955 


1 


4 


1914 




0.572 


0.947 


1 


4 


1918 




0.572 


0.922 


1 


4 


1968 




0.662 


0.968 


1* 


4 


1969 




0.642 


0.965 


1 


4 


1972 




0.572 


0.969 


1 


4 


1973 




0.600 


0.966 


1 


4 


1974 




0.572 


0.962 


1 


4 


1980 


SELfxs, SEL(X4, X^ , Xo), X\ © Xo) 


0.571 


0.966 


1 


3 


2016 


SEL(xi, X2 A (x3 V X4), a;2 A (a;3 V X4)) 


0.773 


1 


1* 


4 


2017 




0.640 


0.972 


1 


4 


2018 




0.710 


0.977 


1 


4 


2019 




0.667 


0.986 


1 


4 


2022 


SEL(x3, SEL(a;2, x\^ X4), SEL(xi, X2, X4)) 


0.661 


0.965 


1 


3 


2025 




0.571 


0.966 


1 


4 


2032 




0.727 


0.971 


1 


4 


2033 




0.667 


0.980 


1 


4 


2034 


SEL(x2, xi, SEL(x4, X3, x\)) 


0.667 


0.980 


1 


3 


2040 




0.600 


0.946 


1 


4 


5736 


EXACT2 


0.572 


1 


1* 


4 


5737 


SYM(0,0,1,0,1) 


0.500 


0.962 


1 


4 


5738 




0.563 


0.975 


1 


4 


5739 




0.500 


0.980 


1 


4 


5742 




0.572 


0.947 


1 


4 


5758 




0.572 


0.922 


1 


4 


5761 




0.500 


0.860 


1 


4 


5763 




0.500 


0.907 


1 


4 


5766 




0.500 


0.936 


1* 


4 


5767 




0.500 


0.933 


1 


4 


5769 




0.500 


0.907 


1 


4 


5771 




0.500 


0.946 


1 


4 


5774 




0.561 


0.966 


1 


4 


5782 




0.500 


0.962 


1 


4 


5783 




0.500 


0.954 


1 


4 


5784 




0.500 


0.946 


1 


4 


5785 




0.500 


0.933 


1 


4 


5786 




0.500 


0.964 


1 


4 


5787 




0.500 


0.955 


1 


4 


5790 


SEL(x3, SEL(x4, xi, X2), a;2 © X4) 


0.500 


0.980 


1 


3 


5801 




0.500 


0.933 


1 


4 
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ID 



Function 



1 query 



2 queries 



3 queries 



5804 
5805 
5820 
5865 
6014 
6030 
6038 
6040 
6042 
6060 
6120 
6375 
6625 
6627 
6630 
7128 
7140 
7905 
27030 



SYM(0,0,1,1,0) 
SEL(x3, SEL(x4, xi,X2), SEL(3;4, X2,xi)) 



SEL(x4, SEL(x3, xi, X2), SEL(3;2, X3, xi)) 

xi © MAJ(a;2, 2:3, 2:4) 
xi © -iNAE(x2, X3, X4) 



Sorted input bits [1] 
xi © SEL(x4, a;2, a;3) 

PARITY 



0.545 
0.500 
0.529 
0.500 
0.600 
0.667 
0.500 
0.572 
0.572 
0.600 
0.667 
0.500 
0.500 
0.500 
0.500 
0.854 
0.500 
0.500 
0.500 



0.961 
0.955 
0.955 
0.954 
0.874 
1 

0.980 
0.951 
0.962 
0.966 
1 

0.900 
0.946 
0.955 
0.954 

1 

1 

0.900 

1* 



1 
1 
1 
1 
1 

1* 
1 
1 
1 
1 

1* 
1* 
1 
1 
1 

1* 
1* 
1* 
1* 



A. 2 Symmetric functions of 5 bits 



Function 


1 query 


2 queries 


3 queries 


4 queries 


(0,0,0,0,0,1) 


0.693 


0.925 


0.988 


0.999 


(0,0,0,0,1,0) 


0.516 


0.761 


0.972 


1 


(0,0,0,0,1,1) 


0.640 


0.798 


0.974 


1 


(0,0,0,1,0,0) 


0.530 


0.616 


1 


1 


(0,0,0,1,0,1) 


0.500 


0.593 


0.995 


1 


(0,0,0,1,1,0) 


0.546 


0.758 


1 


1 


(0,0,0,1,1,1) 


0.600 


0.728 


1 


1 


(0,0,1,0,0,1) 


0.500 


0.640 


0.988 


1 


(0,0,1,0,1,0) 


0.500 


0.517 


1 


1 


(0,0,1,0,1,1) 


0.500 


0.534 


1 


1 


(0,0,1,1,0,0) 


0.600 


0.874 


1 


1* 


(0,0,1,1,0,1) 


0.500 


0.856 


0.998 


1 


(0,0,1,1,1,0) 


0.616 


0.762 


0.969 


1 


(0,1,0,0,0,1) 


0.500 


0.728 


0.967 


1 


(0,1,0,0,1,0) 


0.500 


0.860 


1 


1* 


(0,1,0,1,0,1) 


0.500 


0.500 


1* 


1* 


(0,1,0,1,1,0) 


0.500 


0.616 


0.998 


1 


(0,1,1,0,0,1) 


0.500 


0.784 


0.998 


1 


(0,1,1,1,1,0) 


0.736 


0.962 


0.996 


1* 
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A. 3 Symmetric functions of 6 bits 



Function 


1 query 


2 queries 


3 queries 


4 queries 


5 queries 


(0,0,0,0,0,0,1) 


0.663 


0.900 


0.980 


0.997 


0.9999 


(0,0,0,0,0,1,0) 


0.511 


0.684 


0.940 


0.993 


1 


(0,0,0,0,0,1,1) 


0.640 


0.738 


0.946 


0.993 


1 


(0,0,0,0,1,0,0) 


0.516 


0.572 


0.878 


1 


1 


(0,0,0,0,1,0,1) 


0.500 


0.541 


0.875 


0.999 


1 


(0,0,0,0,1,1,0) 


0.527 


0.751 


0.904 


1 


1 


(0,0,0,0,1,1,1) 


0.589 


0.710 


0.901 


1 


1 


(0,0,0,1,0,0,0) 


0.530 


0.616 


1 


1 


1* 


(0,0,0,1,0,0,1) 


0.500 


0.614 


0.980 


0.997 


1 


(0,0,0,1,0,1,0) 


0.500 


0.504 


0.946 


1 


1 


(0,0,0,1,0,1,1) 


0.500 


0.525 


0.952 


1 


1 


(0,0,0,1,1,0,0) 


0.546 


0.667 


0.864 


1 


1 


(0,0,0,1,1,0,1) 


0.500 


0.625 


0.860 


1 


1 


(0,0,0,1,1,1,0) 


0.556 


0.721 


0.905 


1 


1 


(0,0,1,0,0,0,1) 


0.500 


0.583 


0.882 


0.997 


1 


(0,0,1,0,0,1,0) 


0.500 


0.527 


0.839 


1 


1 


(0,0,1,0,0,1,0) 


0.500 


0.541 


0.843 


1 


1 


(0,0,1,0,1,0,0) 


0.500 


0.517 


1 


1 


1* 


(0,0,1,0,1,0,1) 


0.500 


0.500 


0.985 


1 


1 


(0,0,1,0,1,1,0) 


0.500 


0.520 


0.940 


1 


1 


(0,0,1,1,0,0,1) 


0.500 


0.712 


0.867 


1 


1 


(0,0,1,1,0,1,0) 


0.500 


0.513 


0.840 


1 


1 


(0,0,1,1,1,0,0) 


0.616 


0.762 


0.969 


1 


1* 


(0,0,1,1,1,0,1) 


0.500 


0.722 


0.965 


1 


1 


(0,0,1,1,1,1,0) 


0.625 


0.702 


0.939 


0.992 


1 


(0,1,0,0,0,0,1) 


0.500 


0.652 


0.934 


0.992 


1 


(0,1,0,0,0,1,0) 


0.500 


0.728 


0.967 


1 


1* 


(0,1,0,0,1,0,1) 


0.500 


0.500 


0.836 


1 


1 


(0,1,0,0,1,1,0) 


0.500 


0.667 


0.863 


1 


1 


(0,1,0,1,0,0,1) 


0.500 


0.500 


0.934 


1 


1 


(0,1,0,1,0,1,0) 


0.500 


0.500 




1* 


1* 


(0,1,0,1,1,1,0) 


0.500 


0.553 


0.880 


0.999 


1 


(0,1,1,0,0,0,1) 


0.500 


0.758 


0.908 


1 


1 


(0,1,1,0,1,1,0) 


0.500 


0.640 


0.988 


1 


1* 


(0,1,1,1,1,1,0) 


0.693 


0.925 


0.988 


0.999 


1* 



B Source code 

The following is an example of how the CVX package [11] can be used to determine quantum query 
complexity In this case, we calculate the minimal error probability over all quantum algorithms 
using 2 queries to compute some function / : {0, 1}^ ^ {0, 1} ( given as a column vector). 

cvx_begin 
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cvx_precision best; 



% variables in_i~j : <= i <= n, <= j <= t-1 

variable m00(8,8) symmetric; variable ml0(8,8) symmetric; 

variable m20(8,8) symmetric; variable m30(8,8) symmetric; 

variable ni01(8,8) symmetric; variable ml 1(8, 8) symmetric; 

variable m21(8,8) symmetric; variable m31(8,8) symmetric; 

variable g0(8,8) symmetric; variable gl(8,8) symmetric; 

variable epss; 

minimise ( epss ); 

subject to 

y. Input condition. 

mOO + mlO + m20 + m30 == ones (8, 8); 

y. Running conditions (between 1 and t-1) . 

mOl + mil + m21 + mSl == EO .* mOO + El .* mlO + E2 .* m20 + E3 .* m30; 
y. Output matches last but one query. 

gO + gl == EO .* mOl + El .* mil + E2 .* m21 + E3 .* m31; 

% Output constraints. 
diag(gO) >= (l-epss)*(l-f) ; 
diag(gl) >= (l-epss)*f; 

7, Semidefinite constraints. 

mOO == semidef inite(8) ; mlO == semidefinite (8) ; 
m20 == semidef inite (8) ; m30 == semidef inite (8) ; 
mOl == semidef inite (8) ; mil == semidef inite (8) ; 
m21 == semidef inite (8) ; m31 == semidef inite (8) ; 

gO == semidef inite (8) ; gl == semidef inite (8) ; 

cvx_end 
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